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Abstract 

Phonon transmission across epitaxial interfaces is studied within the lattice dynamic approach. 
The transmission shows weak dependence on frequency for the lattice wave with a fixed angle of 
incidence. The dependence on azimuth angle is found to be related to the symmetry of the boundary 
interface. The transmission varies smoothly with the change of the incident angle. A critical angle 
of incidence exists when the phonon is incident from the side with large group velocities to the side 
with low ones. No significant mode conversion is observed among different acoustic wave branches 
at the interface, except when the incident angle is near the critical value. Our theoretical result 
of the Kapitza conductance Gk across the Si-Ge (100) interface at temperature T = 200 K is 
4.6 x 10 8 WK -1 m~ 2 . A scaling law Gk oc T 2 ' 87 at low temperature is also reported. Based on the 
features of transmission obtained within lattice dynamic approach, we propose a simplified formula 
for thermal conductance across the epitaxial interface. A reasonable consistency is found between 
the calculated values and the experimentally measured ones. 

PACS numbers: 66.70.+f, 44.10.+i 
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I. INTRODUCTION 



The interfacial thermal conductance plays a critical role in nanometer-scaled devices^. 
Kapitza resistance^ was first discovered in 1940s. Although much work, both theoretical and 
experimental,— "2i&2jL&L§j2i has been done since then, the characteristic behavior of phonon 
transmission across the interface is still not clear. The continuum elastic wave model 4 
may be inaccurate when the details of atomistic structures and phonon dispersions are 
considered. The diffusive scattering theory^ neglects the wave property of phonons. Here 
we study phonon transmission across interfaces with the lattice dynamic approach^, which 
can simulate the phonon transmission across interfaces atomistically from the first-principles. 
We concentrate on the phonon transmission at a single epitaxial interface to capture salient 
features of transmission. 

The paper is organized as follows. We first present the method employed in this study. 
Using this method we calculate the dependence of phonon transmission on frequency, the 
azimuth angle, the incident angle, the mode conversion at the Si-Ge interface and at the 
Si-GaP interface in Section II. Some features of phonon transmission are obtained and 
discussed. Finally, on the basis of these features, we propose a simplified formula for thermal 
conductance across interface. A comparison between the calculated thermal conductance 
and the experimentally measured values is made. 



II. METHOD 



We consider the problem of phonon transmission across an epitaxial interface. Two types 
of important crystal structures in semiconductors, diamond and zinc-blend, are chosen in 
this study. The crystals on each side of the interface are assumed semi-infinite^^. For 
simplicity, we also neglect the defects possibly existing between the two crystal solids and 
assume that the interface is an ideal epitaxial interface. Lattice dynamic approach also 
requires that the harmonic approximation for solids is valid. This condition is usually 
satisfied at low temperatures. For the scattering of waves at interface, we remark that 
the nonlinear nonlinear effect may be less important than the elastic scattering at low 
temperature because the thickness of boundary area at the interface is smaller than the 
mean free path of phonons. 
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We follow the scattering boundary method in Ref. QQ. We write the solutions on the 
incident side and on the transmitted side with the unknown component coefficients expressed 
as follows. If a normal mode uf in (u, q) is incident from the left lead, the scattering solution 
for the perfect leads can be assumed as 

n' 

u fi = / i tn"J*U,ri' ( U > In" ) > (lb) 

n" 

where z, Z denote the zth atom in the I unit cell, and n, n' and n" refer to the different 
polarized branches of incident, reflected and transmitted waves. Frequency is denoted as u. 
Wavevectors for the incident, reflected, and transmitted waves are q, q', and q", respectively. 
The superscript L and R indicate the left and the right. In these equations, tffi n , tffi n are 
the amplitude transmission/reflection coefficients from mode n on the lead L to mode n" on 
lead R, and to mode n' on lead L. The wave vectors q' and q" satisfy uj = ov(q') = w n «(q"). 
Note that frequency does not change because the system is linear. 

In this paper, we consider the phonon transmission across two types of epitaxial interfaces: 
the interface between silicon and germanium and the interface between silicon and gallium 
phosphide. Now the problem is to determine the wave numbers for the possible branches 
of the reflected and the transmitted waves. Since the system is homogeneous in the x 
and y directions, the transverse components of the wave vector for the reflected waves q' 
and for the transmitted waves q" have the same values as that of the incident wave q, 
that is q' x = q" = q x ; q' y = q'y = q y . The longitudinal components, q' z and q", satisfy 
u>(q) = w(q') = co>(q"). We solve the nonlinear equation to find q' z or q" for the (111) interface 



with the numerical method in Ref. 
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j. For the (100) interface, we utilize a more efficient 
eigenvalue method to obtain the solutions of nonlinear equations.- The wave vectors q' z and q" 
are identified from these solutions to satisfy the following conditions: (i) The reflected waves 
should have the negative group velocity f 2 (q') so that the phonon energy of the reflected 
waves will propagate back to the — z direction, (ii) The transmitted wave should have the 
positive group velocity so the phonon energy of the transmitted waves will propagate to 
the +z direction. (iii)When the wave numbers for the reflected and transmitted waves are 
complex, they are identified from the solutions to satisfy Im(q' z ) < and Im(q z ) > so 
that these waves decay in space. Actually these decay waves with complex wave number 
do not propagate any energy and the phase difference along z direction for atoms in the 



neighboring unit cells is ir. But they are physically possible modes at the interface and must 
be considered in the scattering boundary method. 




(100) 




(111) 



FIG. 1: Atomic interfacial structure of the (100) and (111) interface. 



Now we consider the dynamic equations for atoms at the interface. The (100) and (111) 
interface between the diamond structure or the zinc-blend structure is illustrated in Fig. [U 
We assume that the transverse directions at the interface are infinitely large so that the 
boundary atoms on each side, e.g., A atoms in the figure, are translational invariant. The 
dynamic equations for A position atoms are equivalent, and similarly for B position atoms 
on the other side of the interface. So we only need to consider the dynamic equations for 
atom A and for atom B at the interface. Under the harmonic approximation, these dynamic 
equations at a given frequency uj can be written as — mjC<j 2 Uj + Yl Kjj ' u j = 0, where 
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mi denotes the atomic mass for each atom, Uj is the oscillation amplitude and Kjj is the 
force constants. The force constants across the interface are chosen as the symmetrized 
Ky = {Kh + Klf)/2. After substituting Eq.([T]) into this dynamic equation, we get six 
equations with six unknown coefficients. This system of equations can be solved by a 
conventional method. The energy transmission from mode (L, n) to mode (R, n") is given 
by 

D R 

(2) 



rpRL 

n"n 



\,RL \2 V n" 
IV'nl ~L 



Here v% is the reduced group velocity^ along z direction for mode n in the left lead, defined 
by = V n/^L where is the group velocity and II is the lattice constant for the left. 



Similar meaning holds for . The total reflection and transmission coefficients 7^ L for 
(L, n) are given, respectively, by 



'^n — / , I Vnl ~ - L > J n ~/ , IVnl ~ L ' 

-R& + Ti = \. (3b) 
The group velocity along 2 direction is calculated through the dynamic matrix as 



(4) 



1 e 1 j=-e 

L aq z 



2 u> eT ■ e 

where D is the dynamic matrix and e the eigenvector of the dynamic matrix-. Note that 
the energy conservation relation, Eq. (13bj) . is satisfied automatically. This can be used as a 
check. With the relation Eq. we can get the Kapitza conductance Gk as 

Gk = ^^^„(q)<(q)T„(q,^„)^^^, (5) 

q,n 

where 7^(q, u n ) is the transmission 8 calculated from Eq. (J3j), is the volume and f(u n , T) is 
the Bose-Einstein distribution for the n — th branch mode. When using Eq. (j3j), we compute 
Kapitza conductance from the Si side. The lattice constants for Si conventional unit cell is 
a = 5.43 A. 



III. MODELS AND NUMERICAL RESULTS 



In this section, we report the dependence of phonon transmission on frequency, the az- 
imuth angle, and the incident angle across the Si and Ge interface. We then consider the 
mode conversion problem at interface. Phonon transmission across a kind of zinc-blend 
interface: Si-GaP is also studied to investigate the lattice structure effect on transmission. 
The atomic masses for Si and Ge are 28amu and 72.61 amu, respectively. After optimizing 
the structure using Tersoff potential, we obtain the linearized force constants under small 
displacement. Due to Tersoff potential truncation function,— only four nearest atoms need 
to be considered for each atom. The phonon dispersions for the Si and Ge along the V — L 
direction can be calculated through the linearized force constants from the dynamic matrix 
and are illustrated in Fig. [2j The maximum frequency of LA branch along V — L for Si 
and Ge are 343 cm -1 and 196 cm -1 , respectively. Compared with the experimental Si and 
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FIG. 2: (A) Phonon dispersion for Si and Ge along T — L and (100) direction. (B) The energy 
transmission T n (to) as a function of frequency uj across the Si-Ge (100) interface for LA waves at 
different angles of incident from Si to Ge. 

Ge phonon dispersions^, there is a discrepancy between the calculated phonon frequency in 
Fig. [2] and the experimental values. This inaccuracy comes from the lack of more neighboring 
atoms for Tersoff potential, apart from the first neighboring ones. The ab initio method 13 
which can include higher order neighboring atoms is desirable for the more elaborate phonon 
dispersion relations. However, we find that the inclusion of higher order neighboring atoms 
will make the calculation of energy transmission very difficult because the solutions of q' z 
or q" at a given frequency u is rather complicated. We estimate that such errors of phonon 
dispersion relation in the Briullouin boundary will not have much influence on the energy 
transmission across the interface. The force constants from the first neighboring atoms play 
a main role in the computation of phonon dispersion relations. The cut-off frequency of 
energy transmission across the interface will not approach closely to the boundaries of the 
Briullouin boundary apart from a few high symmetry directions. The transmission decreases 
very rapidly when the frequency is near cut-off, which can be seen from Fig.2(B). Therefore, 
we think that the inaccuracy of the Briullouin boundary values from Tersoff potential will 
not significantly change the results of energy transmission across the interface. 
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A. Dependence on angular frequency 

The transmission's dependence on angular frequency across a (100) Si-Ge interface is 
shown in Fig. [21(B). The incident acoustic wave is longitudinally polarized with an incident 
angle 9 in xz plane. For waves with the incident angle 9 below about 63°, the transmission 
falls into a very narrow range approximately from 0.91 to 0.98 and decreases slightly with 
the increase of frequency before it reaches its cut-off frequency. Here the cut-off frequency 
means that beyond this point the phonon transmission equals zero. There is an abrupt 
decrease in phonon transmission when the frequency is near this value. When the incident 
wave is normal to the interface along (100) direction from Si to Ge, the transmitted waves 
are also along (100) direction. Thus, the cut-off frequency of transmission can be easily 
determined from the phonon dispersion of Ge. This can be shown from Fig. [2]^B). The 
cut-off frequency is about 197 cm -1 for the wave with the incident angle 9 = 0, while the 
endpoint frequency of the longitudinal acoustic wave branch for Ge along (100) direction is 
200 cm -1 , which agrees with the cut-off frequency of transmission for the normal incident 
wave. In Fig. [2](B), for clarity, we only plot the transmission before its abrupt jump to 
zero. For transversely polarized incident acoustic waves, we also have observed a similar 
phenomenon. It can be concluded that the acoustic phonon transmission from Si to Ge 
varies in a narrow range with the increase of frequency, independent of the the property 
of its polarization. Note that the phonon wave is incident from the material (Si) with a 
larger group velocity to the material (Ge) with a smaller group velocity. This conclusion 
also holds in our calculations of the transmissions for other incident acoustic waves, not only 
in xz plane, but also with azimuth angle (p. We remark that this small change of phonon 
transmission with the frequency for incident wave at the incident angle 9 can be neglected 
for a approximate estimation of thermal conductance across an interface. 

B. Spatial angular dependence 

1. Azimuth angle and symmetry of the interface 

We find that the phonon transmission's dependence on the azimuth angle of the incident 
wave is related to the symmetry of the interface. The results of the transmissions for the 
longitudinal acoustic waves incident from Si to Ge are shown in Fig. [31 It can be seen 
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from Fig. [3] that the phonon transmission shows little dependence on the azimuth angle 
(p at incident angle 9 < 45° both for the (100) and (111) interface. However, when the 
incident angle increases beyond 45°, there are peaks in transmission with the variation of 
the azimuth angle. The peaks reveal the anisotropy of the interface. There is a four-fold 
symmetry axis in the [100] direction in Si or Ge and a three-fold symmetry axis along the 
[111] direction. These symmetries have evidently been shown in the phonon transmission 
across the related interface. In Fig. [3], there are four- fold symmetrical peaks separated by 7r/2 
in the transmission across the (100) interface. Three three-fold symmetrical peaks separated 
by 27r/3 appear in the transmission across the (111) interface. 
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FIG. 3: Dependence of the longitudinal acoustic phonon transmission on the azimuth angle if. 
All transmissions plotted are at angular frequency u) = 15 cm -1 and are incident from Si to Ge. 
Different lines in the figure correspond to the transmissions for wave with different incident angle 
9 as the indicated in the figure. Left figure shows the results for the (100) interface, right figure 
for the (100) interface. 



2. Critical incident angle 

We next report the result of the dependence of energy transmission on the incident angle. 
A continuum wave incident on a surface with an incident angle 9 is refracted in accordance 
with Snell's law. It would be an interesting question if there is a critical angle for discrete 
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lattice waves. We calculated the phonon transmission incident from Si to Ge and incident 
from Ge to Si. The dependence of transmission on the incident angle are illustrated in Fig. H] 
for the (100) Si-Ge interface. It can be seen from Fig.H]that when the angle 9 for the phonon 




FIG. 4: Dependence of phonon transmission on the incident angle 6. The left figure shows the 
results of longitudinal acoustic waves (LA) and the right figure is for the transverse polarized 
waves (TA). The solid lines indicate the transmission from Si to Ge and the dashed lines denote 
the transmission from Ge to Si. All the transmissions by lattice dynamic approach are calculated 
at angular frequency u = 15 cm" 1 . The dotted line in the left figure is the result calculated from 
the continuum wave model. 

incident from Si to Ge increases, energy transmission decreases slowly. The incident angle 
from Si to Ge can be extended to as large as 90°. In contrast to the transmission from Si 
to Ge, for waves incident from Ge to Si, there exists a critical angle, about 38° for LA and 
38.5° for TA, above which the transmission is very small. We can estimate the critical angle 
with the help of Snell's law for continuum wave. The group velocities calculated from the 
dynamic matrix in this paper for longitudinally polarized waves along the [001] direction 
are t>£* ~ 6.87 Km/s and v1 e ~ 3.78 Km/s, repectively. The Snell's law of continuum wave 
model gives the critical angle 9 C = sin -1 (f^ e /t>£*) = 33.4° from Ge to Si. This value is 
a little lower than the observed critical angle value. However, we can still think that the 
Snell's law holds approximately. 

When the incident waves from Si to Ge and from Ge to Si are both normal to the 
surface (6 = 0), the transmissions from both sides are 0.98 for LA waves and 0.95 for 
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TA waves. We use the acoustic mismatch model^ to estimate AZsiZcel '{Z$i + Zee) 2 ~ 
0.97 for LA, 0.94 for TA, where Zs% and Zc e denote the acoustic impedance defined as Z = 
pv. Here p is the mass density and v is the group velocity. Here we take p$i = 2.329 x 
10 3 Km/s and pQ e = 5.323 x 10 3 Km/s. The values of group velocities for LA waves are 
stated in the previous paragraph. It can be seen that the acoustic mismatch model in Ref |4 
well describes the transmission for normal incident waves. 

To investigate the difference in the transmission on the incident angle between the lattice 
dynamic approach and the continuum wave model, we also calculated the energy transmis- 
sion of the continuum wave.— The results for Si-Ge interface is illustrated in Fig. HI The 
procedure of calculation is as follows. First, the amplitudes transmission^ is computed from 
the Fresnel equation t = ZCO sg^l os g> ? where z = y/pc and z' = y/p'd is the wave impedance 
for the incident and refraction crystal. Here the mass density p and p' take the value in the 
preceding graph. The stiffness constants c, d are 0.796, 0.680 x 10 11 N/m , respectively. The 
angle of incidence 9 and the angle of transmission 9' satisfy the Snell's Law. The energy 
transmission for the continuum wave is calculated through the formula T = \t\ 2 p - cos f , 

° I I pv cos a ' 

where v = 8.43Km/s, v' = 4.87Km/s are the group velocities for the incident and refraction 
continuum waves. It can be seen that the isotropic continuum wave gives a similar depen- 
dence of transmission of the incident angle. The result for the continuum wave is close to 
the result of lattices wave. 

C. Mode conversion at the interface 

The mode conversion at the interface would be an interesting problem to pursue. There 
are two types of mode conversion for thermal transmission across the boundary: acoustic- 
optical (AO) and acoustic-acoustic (AA). We neglect the optical-optical conversion because 
the corresponding frequencies do not overlap in our model. We think that this optical-optical 
conversion's contribution to thermal transport is trivial because of its relatively low group 
velocity and high energy. 

AO conversion. We did not observe significant AO conversion in our simulation of phonon 
transmission at the Si-Ge interface. The AO conversions in our results are very small, no 
more than 10% and fall in a very narrow frequency range. This result is different from 



the reported results in Ref. We think this maybe come from the empirical potential 
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chosen in our paper, which results in different group velocities. To estimate the maximum 
ratio of a possible acoustic-optical mode conversion, A simplified ID toy chain model 8 can 
be composed because it can be assumed that the AO conversion will be easier for both 
longitudinal waves in one dimension. However, the results in Ref. ^ find that the energy 
transmission contributed by the AO conversion is trivial in comparison with the acoustic- 
acoustic transmission. This is due to the large mismatch in their group velocities. 
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FIG. 5: Demonstration of mode conversion LA — * LA + TA near the critical angle from Ge to 
Si. (A) The ratio of transmission converted to LA; (B) The ratio of transmission to TA; (C) The 
energy transmission T n (u) as a function of angular frequency uj across the Ge-Si (100) interface 
for LA waves at different angles of incident from Ge-Si. 



AA conversion. Apart from the possible AO conversion, there are several different po- 
larized acoustic branches. Can the the conversion among these different acoustic branches 
occur at the interface? Our simulation did not find significant conversion among these dif- 
ferent polarized branches, except when the angle of incidence is near the critical angle. The 
reflected and transmitted waves are both LA waves when the wave incident is longitudinally 
polarized. The same is true for the TA wave. The AA conversion among different polarized 
acoustic branches only takes place when the angle of incident is close to the critical angle for 
wave incident from Ge to Si. Fig. [5] shows the transmission of the longitudinally polarized 
wave incident from Ge to Si. The mode conversion of LA — > LA + TA is shown in pictures 
(A) and (B) in Fig. [5j It can seen that the ratio converted from LA to LA decreases with 
the increase of the angle of incident. The ratio converted from LA to TA increases rapidly 
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to about 0.12 near the critical angle. This behavior of mode conversion taking place near 
the critical angle maybe result from the suppression of transmission due to the Snell law. 
We also find that the presence of the possible reflected waves influence the transmission. 
For example, in contrast to the transmission from Si to Ge, the dependence of LA mode 
wave transmission on frequency from Ge to Si shows a rich character as illustrated in Fig. 
For wave with the incident angle 9 ^ 0, the transmission first decreases with the increase of 
frequency. But when the frequency goes over a certain value, for example to = 93 cm -1 for 
9 = 36°, the transmission begins to increase. This behavior can be understood by the mode 
conversion at the interface. It can be seen from Fig. [2] that the maximum frequency for TA 
modes in Ge is about 98 cm -1 . When the frequency is below this value, there are TA modes 
for the reflected wave; but over this value, the reflected wave cannot be converted into TA 
modes. The transmission increases due to lack of reflected modes, . 

D. Transmission across the interface between diamond and zinc-blend structure 

During the calculation of the transmission of Si-Ge interface, the atomic masses are 
uniformly distributed both for Si and Ge on each side. There is another important type 
of lattice structure in the semiconductor materials. That is the zinc-blend structure, where 
the atomic masses in the sublattice are different. Here we consider an example of gallium 



1.0- 




FIG. 6: Left figure: Phonon dispersion for acoustic waves in GaP. Right figure: Dependence of 
phonon transmission for LA waves on the incident angle 9 at angular frequency u = 15 cm -1 from 
Si to GaP. 
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phosphide(GaP). To simplify the problem, we only consider the valence-force and neglect the 
Coulomb interaction between the charges when computing the phonon dispersion of GaP. 
This is a valid approximation when we consider the acoustic phonons that play major role 
in the phonon transmission across the interface. The calculated acoustic phonon dispersion 
curves along V — L for Gap is plotted in Fig. [6J The endpoint values of frequencies along 
r — L for LA and TA are about 208 cm -1 and 116 cm -1 , respectively. It can be seen that 
valence-force constants approximately describe the acoustic branches^ of GaP. The results of 
phonon transmission for LA and TA modes across a Si-GaP interface are illustrated in Fig. [61 
We find a consistent dependence on frequency, spatial angular and a similar phenomena of 
mode conversion with that of Si-Ge interface. The dependence of transmission on the angle 
of incidence is illustrated in Fig. [6J 

E. Temperature dependence of Kapitza conductance 

The Kapitza conductance with the change of temperature is calculated using Eq. (GO), and 
are illustrated in Fig. [3(A). The Kapitza conductance for Si-Ge [100] interface calculated 
from our model is Gk = 4.6 x 10 8 WK _1 m~ 2 when T = 200 K. When the temperature goes 
beyond 200 K, we find that the Kapitza conductance changes little with the temperature 
and is saturated. For comparison, we plotted the heat capacity of Si in Fig. [7(B). It can be 
seen that the heat capacity continues to increase with the temperature when T > 200 K. In 
comparison with the heat capacity, the saturation of Kapitza conductance can be accounted 
by negligible contribution of energy transmission from high frequency at low temperatures, 
as illustrated in inset of Fig. [71 The Kapitza conductance scales as T 2,87 , while the heat 
capacity scales as T 3 in accordance with the Debye model. We have sampled enough points 
in the first Brillouin zone to ensure that Kapitza conductance and heat capacity converge 
numerically. However, due to the small deviation from the value of three^we cannot rule out 
the possibility that the exponent for Kapitza conduction is also 3. Ref. |7| reported Kapitza 
conductance scaled T 3 at low temperature for fee interface irrespective of the properties for 
the left and right lead. The temperature dependence of Kapitza conductance is an intriguing 
problem 2 -, though much experimental work has been done on this field. Most experiments 
gave T a with a < 3 for solid interface as reviewed in Ref. |2j. So far no experiment result is 
available for the temperature dependence for the Si-Ge interface. Compared with the results 
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FIG. 7: (A) The temperature dependence of Kapitza conductance for Si-Ge(lOO) interface. (B) 
The corresponding heat capacity of Si. 
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of Ref . |7| our discrepancy from T 3 comes from the anisotropy of the energy transmission 



because of the diamond structure used for calculation of the transmission, while Ref. 
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took 



an isotropic assumption for their calculation. 



IV. A NEW SIMPLIFIED APPROXIMATE FORMULA 

To get a simplified expression, we make the follow approximation on the basis of the 
features of transmission discussed above. (1) For transmission with a given incident angle 
9, we assume that it is independent of the frequency. (2) We neglect the dependence on 
the azimuth angle. (3) The conductance can be calculated from either the left or the right 
side.- They give the same result across the interface. It is found that the transmission from 
the high group velocity to the low velocity side has a simpler feature. We calculate the 
transmission from the high to the low velocity side. The dependence of transmission on 
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the incident angle for a polarized waves is formulate as T n = T^ 1 cos 6. Here 7^™ is the 
transmission at uj = with the angle of incident 8 = 0, and is given by 

T n _ 4Z 1 Z 2 (6 x 

Z? = Pi v -f-, with i = 1,2. (6b) 

Here is the acoustic impedance for the polarized wave and U is the lattice constant. We 
have incorporated the effect of difference in lattice constants of the crystal. It can be seen 
from Fig. H] that the cosine function is a valid approximation. (4) Mode conversion at the 
interface can be neglected. (5) Since thermal conductance is contributed mainly by acoustic 
waves, we can assume that the upper limit of frequency of acoustic waves is characterized by 
the Debye temperature of the side with the group velocity. A linear dependence of phonon 
dispersion relation u = qv n is used because thermal conductance is mainly due to acoustic 
waves. The group velocity v z n along z in Eq. \5\ is given by v z n = v n cos 9. Under these 
conditions, we have a simplified thermal conductance (Kapitza conductance) as 



XD x 4 e x 



Here v n is the group velocity of the incident side with the large group velocity. The upper 
limit xd is given as id = -jr, where 9d is the Debye temperature for the other side with low 
group velocity. We use Eq. (j7j) to calculate a few interfaces across which the experimentally 

TABLE I: Values for mass density p, lattice constant I, the longitudinal group velocity vl, the 
transversal group velocities vt, heat capacity C v and the Debye temperature Op. These values are 



from Ref. 
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p(10 3 Kg/m 3 ) 


/(A) 


vl( Km/s) 


vt{ Km/s) 


a(10 6 J/m 2 K) 


MK) 


Si 


2.329 


5.43 


8.43 


5.84 


1.63 


640 


Al 


2.699 


4.05 


6.24 


3.04 


2.42 


394 


A1 2 3 


3.97 


4.76 


10.89 


6.45 


2.656 


1024 


Bi 


9.79 


4.75 


1.972 


1.074 


1.194 


120 


Pb 


11.598 


4.95 


2.35 


0.97 


1.463 


105 


Diamond 


3.515 


3.57 


17.52 


12.82 


1.8278 


1860 



measured thermal conductance is available. The parameters taken are shown in Table [B 
The results of calculation are shown in Fig. [HJ 
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FIG. 8: The temperature dependence of Kapitza conductance calculated from Eq. [7J 

TABLE II: Comparison of thermal conductance across interface for the experimental measured 
values and the calculated values at T = 200-KT. The experimental values are from Ref. (J The 
values are in the units of 10 8 W/m 2 K. 







A1/A1 2 3 


Pb/Diamond 


Bi/Diamond 


Pb/Si 


Bi/Si 


Exp. 


Values 


1.2 


0.125 


0.06 


0.13 


0.1025 


Cal. 


Values 


2.5 


0.013 


0.02 


0.10 


0.0640 



We list the experimental values and our calculated values in Table ILT1 It can be seen that 
the calculated values are consistent with the experimental values except for the interface 
Pb/Diamond. The large deviation for Pb/Diamond interface between the calculated and 
experimental values maybe come from the nonlinear scattering effect at the interface. It 
cannot be explained in the reach of the method employed in this paper. 



V. DISCUSSION AND CONCLUSION 



We have studied the features of phonon transmission across the epitaxial interfaces by 
lattice dynamic approach. The transmission is found to change slightly with the frequency 
for polarized waves with a given incident angle from the high group velocity side to the 
low group velocity side. The dependence of transmission on the azimuth angle is related 
to the symmetrical properties of the interface for waves with larger angle of incident. A 
critical angle exists for transmissions from the low group velocity side to the higher group 
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velocity side. Dependence of transmission on the mode conversion at the interface is trivial 
except when the incident angle is close to the critical angle. Thermal conductance across the 
epitaxial interfaces is dominantly contributed by the acoustic waves and is insensitive to the 
mass distribution of different lattice positions. Kapitza conductance across the Si-Ge(lOO) 
interface shows a T 2 ' 87 dependence on temperature, which is very close to T 3 behavior. A 
simplified formula for the estimation of thermal conductance across the interface is proposed 
in the light of features found by lattice dynamic approach. We remark that this formula 
can give a valid estimation of Kapitza conductance across solid epitaxial interfaces when 
nonlinear phonon scattering is unimportant. We think nonlinear scattering sometime maybe 
important for the interfaces between materials with the higher Debye temperature and the 
lower Debye temperature, such as Pb-diamond. Such a nonlinear effect cannot be captured 
in the current lattice dynamic approach. From the calculation, the nonlinear effect maybe 
contribute to the increase of thermal conductance across the solid interface by breaking the 
selecting rules between the linear phonon dispersion relations. The phase relations among 
the incident wave and the reflected, transmitted waves have not been considered because 
they will not influence the energy transmission across the interface. A possible existence of 
lattice dislocation or disorder is not included in the present lattice dynamic approach. 
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